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Simulation  of  Leading  Edges  of  Backscatter 
lonograms  for  Ionospheric  Inference 


1.  IMROIH  CTION 

This  report  presents  the  results  of  ray  computations  made  for  simulating  the 
principal  features  of  swept-frequency  backscatter  radar  ionograms,  with  emphasis 
on  the  leading-edge  traces  appearing  in  such  ionograms.  The  eventual  objective  is 
to  infer,  from  observed  backscatter-radar  ionograms,  important  features  of  the 
ionospheric  electron-density  distribution  in  the  observational  sector  of  the  radar, 
and  thereby  infer  the  major  effects  of  the  ionosphere  upon  the  performance  of  the 
radar.  The  eventual  procedure  is  to  start  with  an  initially  assumed,  trial  electron 
distribution,  and  then  iteratively  alter  the  electron  distribution  until  the  computed 
leading-edge  traces  approximately  match  the  corresponding  traces  in  the  given  set 
of  observed  backscatter  ionograms. 

The  ray  computations  were  made  under  conditions  of  three-dimensional  varia- 
tions in  the  ionosphere  at  polar  latitudes.  Figure  1 shows  a shaded  area  where 
auroral  activity  and  field-aligned  ionization  were  observed  near  local  midnight 
hours  on  9 December  1971.  (The  southern  edge  of  the  shaded  area  constitutes  the 
southern  boundary1  of  the  auroral  oval. ) At  this  time,  the  center  line-of-sight  of 

(Received  for  publication  24  November  1976) 

1.  Buchau,  J. , Pike,  C.  P.,  and  Wong,  M.  S.  (1973)  Detailed  Specification  of  the 
Arcttc  Ionosphere  and  an  Application  to  Three-Dimensional  Ray  Tracking. 
AFCRL-TR-73-0726. 
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the  radar  grazed  the  boundary  between  the  aurorally  active  and  F-layer  trough 
regions.  In  this  situation,  it  is  to  be  expected  that  the  radar  antenna  beam  along 
the  center  line-of-sight  azimuth  would  experience  strong  off-great -circle 
deviations. 


Figure  1.  I.ines  of  Sight  from  Backscatter  Hadar,  and  Location  of  Field-aligned 
Ionization  Observed  in  Shaded  Area  near  2300  LT  on  9 December  1971.  (After 
Buchau  et  al. ) The  map  is  plotted  in  accurate  geomagnetic  latitude  and  longitude 
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In  ray  computations  using  the  Jones  program  for  solving  Hamiltonian  ray 
equations  for  anisotropic  media,  the  dipolar  geomagnetic  coordinate  system  is 
most  advantageously  used.  As  seen  in  Figure  1,  accurate  geomagnetic  coordinates 
are  appropriate  for  describing  important  ionospheric  features  at  polar  latitudes. 
Thus,  it  is  necessary  to  use  the  geographic,  the  dipolar-geomagnetic,  and  the 
accurate-geomagnetic  systems  of  coordinates  interchangeably  in  ray  and  ionogram 
computations. 

Figure  2 shows  the  transformation  between  geographic  and  accurate-geomag- 

3 

netic  coordinates.  The  contours  of  geographic  colatitude  are  shown  in  1-degree 
increments.  Origin  of  the  accurate  geomagnetic  (polar)  coordinates  is  at  the  inter- 
section of  the  row  and  column  of  dots  spaced  at  5-degree  increments  of  geomagnetic 
colatitude.  Figure  3 shows  the  transformation  between  dipolar-geomagnetic  and 
accurate -geomagnetic  coordinates.  The  curves  in  Figure  3 were  obtained  by  using 

4 

a computer  subroutine  developed  by  Bhavnani  for  inverting  the  Gustafsson 
tables,  5 which  converts  geographic  to  accurate-geomagnetic  coordinates,  and  by 
using  a transformation  matrix  between  geographic  and  dipolar-geomagnetic 
coordinates. 

For  convenience,  as  in  a previous  report,  r>  an  arbitrary  electron  distribution 
is  deliberately  taken  to  consist  of  one  part  extending  over  the  entire  geographic 
region  of  interest,  and  a second  part  over  a smaller  area  (for  example,  strong 
E-laver  ionization  over  latitudinal  and  longitudinal  intervals  associated  with  auroral 
activities).  For  each  of  these  two  parts,  the  electron  density  varies  with  height  as 
a continuous  sequence  of  up  to  five  monotonic  segments  in  the  bottom  side  of  the 
ionosphere:  each  segment  is  a sine-squared  function  with  zero  slopes  at  the  ends, 
increasing  or  decreasing  with  altitude.  TheChapman  function  is  used  for  the  top 
side  of  the  ionosphere. 


2.  Jones,  R.M.,  and  Stephenson,  J.J.  (1975)  A Three-Dimensional  Ray  Tracing 

Computer  Program  for  Radio  Waves  in  the  Ionosphere,  OT  Report  75-76, 

II.  S.  Department  of  Commerce,  Office  of  Telecommunications.  (For  sale 
by  Government  Printing  Office.) 

3.  Whalen,  J.A.  (1970)  Auroral  Oval  Plotter  and  Nomograph  for  Corrected 

Geomagnetic  Local  Time,  Latitude,  and" Longitude,  A FCRL-TR -70 -0422. 

4.  Bass,  J.N. , and  Bhavnani,  K.  H.  (1974)  Analysis  and  Programming  for 

Scientific  Research,  AFCRL-TR-74-O4T50. 

5.  Gustafsson,  G.  (1970)  A revised  corrected  geomagnetic  coordinate  system, 

Arkiv  Geofysik  5:595-617. 

6.  Wong,  M.S.  (1975)  Iterative  Ray-Tracing  Simulation  of  Minimum  Group-Path 

Traces  in  Swept -Frequency  Backscatter  Ionograms,  AFCRL-TR-75-0241. 
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Figure  2.  Curves  of  Geographic  Colatitudes  and  Longitudes  Plotted  on  Polar  Co- 
ordinates for  Accurate  Geomagnetic  Latitude  and  Longitude.  (After  Whalen.) 


Figure  4 illustrates  how  the  sine-squared  segments  may  be  used  for  the  two 
parts  of  the  electron  distribution.  For  each  segment,  the  maximum  and  minimum 
electron  densities,  as  well  as  their  altitudes,  all  may  be  functions  of  colatitude 
«ac  ln  the  accuratc  geomagnetic  coordinate  system.  As  shown  in  Figure  5a,  each 
function  of  0ac  consists  of  a continuous  sequence  of  sine-squared  function  segments 
with  zero  slopes  at  the  end  points  of  each  segment. 
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DIPOLAR  COLATITUDE  (DEG) 


Figure  3.  Curves  of  Accurate  Geomagnetic  Latitudes  and  Longitudes  Plotted  on 
Cartesian  Coordinates  for  Dipolar  Geomagnetic  Latitude  and  Longitude 


J I I 1 I L-L  1 

N2  N3  N4  N5  N|  N3 

ELECTRON  DENSITY 


Figure  4.  Two  Separate  Electron-density  vs  Height  Profiles  Lsed  for  Representing 
any  Given  Electron  Distribution.  Each  profile  consists  of  five  sine-squared  seg- 
ments joined  continuously  for  bottom  side  of  ionosphere 
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Figure  5a.  Variation  with  Accurate  Geomagnetic  Colatitude  of  Ci  n- 
ical  Frequencies,  and  Critical-frequency  Heights  of  Sine-square  Seg- 
ments Shown  in  Figure  4 


ACCURATE  MAG  COLATITUDE  (DEG) 

Figure  5b.  An  Alteration  of  the  Critical-frequency 
Function  Fa  in  Figure  4 with  the  Dash-line  Segment 
Replacing  the  Original,  Steeper  Solid-line  Segment 
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Figure  fi.  The  Critical  Frequencies  and  Heights  Fa,  F5,  F^,  II4  and  H5  in  Figure  5 
Plotted  as  Functions  of  Ground  Range  from  the  Backscatter' Radar  in  Various  Azi- 
muths: (a)  20°  and  25°;  (b)  30  and  3 5°;  (c)  40°  and  45°;  (d)  50°  and  55° 


13 


-600 

-400 

H 

0 200“ 
UJ 

1 0„ 


DIPOLAR  AZIMUTH  ■ 400 


1000  2000  3000  4000  5000 


DIPOLAR  AZIMUTH  ■ 45  0 


2000  3000  4000  5000 


14 

5 

2r 
0- 


o 

UJ 

cr 
u_ 

_j 

< 

o 

O ' 
°0 

(c) 


/ 


L 


\ 


V 


1000  2000  3000  4000  5000 


1000  2000  3000  4000  5000 

GROUND  RANGE  (km) 


6 

4 

2 

°o 

2 1 
0j- 


1000  2000  3000  4000  5000 


1000  2000  3000  4000  5000 


GROUND  RANGE  (km) 


-600 

c |_ 

H5  DIPOLAR  AZIMUTH  • 

500 

600  r 

u _ 

DIPOLAR  AZIMUTH  ■ 55  0 

~400i 

y~\ 

400: 

h5 

/~\ 

t L \- 

<5200; 

UJ  I 

1 n 

h4 

200 j 

h4 

°0 

1000  ' 2000 

3000  ' 

4000 

5000 

°0 

1000 

2000 

3000  4000  5000 

6 

6r 

X4 

4- 

5 i 

F5 

l 

F5 

O 2* 

2 r 

\_V 

UJ  • 

h 

<r  n 

^ 0 

1000  " 2000 

3000 

4000 

5000 

°0 

1000 

2000 

3000  4000  5000’ 

< 

u 

E2[ 
O i 

h 

2[ 

F5 

°0~ 

1000  2000 

3000 

4000 

5000 

°o” 

1000 

2000 

3000  4000  5000 

(d) 

GROUND  RANGE  (km) 

GROUND  1 

RANGE  (km) 

Figure  6.  The  Critical  Frequencies  and  Heights  F4,  F5,  F5,  H4  and  in  Figure  5 
Plotted  as  Functions  of  Ground  Kange  from  the  Backscatter' Hadar  in  Various  Azi- 
muths: (a)  20°  and  25°;  (b)  30°  and  35  ; (c)  40°  and  45°:  (d)  50°  and  55°  (Cont.  ) 
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For  the  results  presented  in  the  report,  the  curves  in  Figure  5a  were  taken 
to  remain  unchanging  over  a 90-degree  sector  in  accurate  geomagnetic  longitude 
for  local  times  near  midnight.  However,  the  electron-density  processing  sub- 
routine used  does  include  provisions  for  perturbational  variations  with  <j>  of  the 
•j  3.C 

curves  in  Figure  5a. 

Figure  5b  shows  an  example  of  the  major-feature  alterations  that  can  be  made 

7 

very  expeditiously,  using  an  electron-density  preprocessing  subroutine,  on  the 
function  F_(p)  appearing  in  Figure  5a. 

Viewed  as  functions  of  ground  range  from  the  transmitter  at  ( f)Q,<t>0 ) (32.0, 

1.8)  degrees,  and  in  various  chosen  directions  of  initial  ray  azimuth  reckoned  in 
the  dipolar-geomagnetic  coordinate  system,  the  curves  in  Figure  5a  become  trans- 
formed to  those  in  Figures  6a,  b,  c,  and  d. 


I.  COMM  TKIMOMMiR  Ml  tllUt  M TFKISTICS 


For  the  electron-density  distr'bution  N(r)  represented  by  Figures  4 and  5a, 
ray  trajectories  were  computed  for  various  initial  take-off  angles  in  elevation, 
and  with  frequency  and  initial  azimuth  as  parameters.  Azimuth,  colatitude  f),  and 
longitude  <j>  are  all  reckoned  in  the  dipolar  geomagnetic  system  of  coordinates  in 
this  section. 

Figure  7 shows  the  group-path  lengths  of  rays  from  the  transmitter  to  back- 
scattering  points  either  in  a field-aligned  region  where  the  wave-normal  directions 
intercept  the  accurate  geomagnetic  field  vector  at  90.  0 degrees  in  aspect  angle,  or 
in  a landing  area  of  the  rays  at  ground  level.  The  curves  were  computed  with  the 
transmitter  at  {f)g,  <f>o>  = (32.0,  1.8)  degrees,  for  frequency  = 5.75  MHz,  and  for 
initial  azimuths  of  the  rays  at  10,  15,  ...  55  degrees.  The  larger  (underlined) 
numbers  near  the  slanting  U-shaped  curves  denote  initial  azimuths;  the  smaller 
numbers  denote  initial  elevation  angles  of  the  rays  concerned.  The  horizontal  axis 
represents  the  azimuths  of  back-scattering  points  viewed  from  the  transmitter. 

In  Figure  7,  the  rays  emitted  at  one  azimuth  are  divided  typically  into  sub- 
intervals of  initial  elevation  angle.  Each  subinterval  gives  rise  to  a slanting 
U-shaped  diagram;  more  than  one  U-shaped  diagram  may  belong  to  one  single 
propagation  mode,  such  as  one-hop  E,  one-hop  F,  backscatter  by  field-aligned 
ionization,  or  combination  of  these  modes.  Thus  there  may  be  more  than  one 
group-path  minimum  for  rays  emitted  in  a single  azimuth  and  propagated  via  a sin- 
gle mode.  These  group-path  minima,  plotted  as  functions  of  frequency,  in  turn 
give  rise  to  multiple  "leading-edge"  traces  for  a single  propagation  mode  in 

7.  Langworthy,  B. , and  Barrett,  T.  (1976)  Final  Report  under  Contract 
F 1 9628 -7 6 -C -0029.  RADC/ET,  Hanscom  ATE. 
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Figure  7.  Group-path  Lengths  for  Rays  Transmitted  at  Various  Take-off  Angles 
of  Elevation  and  at  Fixed  Initial  Azimuths,  Plotted  Against  the  Azimuths  at  the 
Backscattering  Points  of  the  Rays  Concerned.  Initial  azimuths  10°,  15°,  ...  55°: 
transmitter  at  (32°,  1.8°);  frequency  5.  75  MHz;  N(r)  shown  in  Figures  4 and  5a. 
Bottom  half  for  ground  backscatter;  top  half  for  backscatter  by  field-aligned 
ionization 


observed  backscatter-radar  ionograms.  For  each  point  on  a U diagram,  the  back- 
scatter azimuth  minus  the  initial  ray  azimuth  measures  the  off-great-circle  path 
deviation  of  the  backscatter  points  concerned  from  the  initial  ray  azimuth;  a nega- 
tive difference  indicates  a clockwise  deviation  in  azimuth.  Off-great-circle  devia- 
tions of  more  than  -20  degrees  occur  for  some  rays.  However,  for  rays  near  the 
group-path  minimum  in  a U diagram  (that  is,  rays  which  have  nearly  maximum 
ray  density),  the  deviation  is  typically  a few  degrees. 

Figures  8,  9,  and  10  show  the  leading-edge  (minimum  group-path)  traces 

obtained  with  the  transmitter  at  (o  , d>  ) (32.0,  1.8),  (34.0,  1.8),  and  (38.0,  1.8) 

o o 

degrees  respectively,  and  for  various  initial  ray  azimuths  between  20  and  50 
degrees.  The  propagation  modes  involved  in  the  curves  shown  are  ground  back- 
scatter via  the  E or  F region  (labeled  EG,  or  FG);  and  backscatter  by  field -aligned 
ionization  in  the  E or  F region  (El, or  FT),  at  points  where  the  aspect  angle  between 
the  wave-normal  direction  of  the  ray,  and  the  direction  of  the  accurate  geomagnetic- 
field  vector  is  90°.  The  curves  drawn  with  dashes  indicate  that  the  initial  elevation 
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F’igurS  9-  Leading- edge  Traces  for  N(r)  Shown  in  Figures  4 and  5a.  Transmitter 
at  (34  , 1.8°);  (a)  initial  az  = 20°;  (b)  Initial  az  3CT;  (c)  initial  az  45°; 
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Figure  10.  Leading-edge  Traces  for  N(r)  Shown  in  Figures  4 and  5a;  (a)  initial 
az  = 20°;  (b)  initial  az  = 30°;  (c)  initiafaz  = 40°;  (d)  initial  az  = 50° 
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angles  of  the  minimum-group-path  rays  are  small  (two  degrees  or  smaller).  Dots 
on  a curve  indicate  that  the  perpendicularity  between  wave-normal  and  geomagnetic- 
field  directions  occurred  south  of  the  auroral  oval  region  where  field -aligned 
ionization  prevails  with  highest  probability.  It  is  seen  that  cusps,  each  occurring 
at  the  intersection  of  two  leading-edge  traces  belonging  either  to  the  same  or  dif- 
ferent propagation  modes,  appear  prominently.  The  leading-edge  traces  shown 
were  determined  by  a computer  program  that  considered  solely  the  minimum  group 
path  length  for  the  slanting  U-shaped  diagrams  involved  such  as  the  ones  shown  in 
Figure  7.  In  practice  these  traces  may  need  to  be  modified  on  consideration  of 
ray  density,  ionospheric  absorption  along  ray  paths,  and  radar  antenna  patterns. 

At  the  group-path  minimum  of  a U diagram,  ray  density  is  usually  close  to  maxi- 
mum with  fair  approximation;  also,  absorption  computed  with  the  Elkins -Rush 

g 

model  for  electron-collision  frequency  in  the  ionosphere  does  not  vary  greatly 
for  neighboring  rays. 

Figures  8,  9,  and  10  show  that  the  principal  features  of  the  leading-edge  traces 
change  gradually  with  increase  of  initial  ray  azimuth  until  the  rays  graze  the  lati- 
tudinal belt  shown  in  Figure  5a  where  the  critical  frequency  F,.  falls  steeply  with 
geomagnetic  colatitude.  Then  the  leading-edge  traces  for  the  E -region  propagation 
modes  disappear,  because  in  Figure  5a,  the  E-region  ionization  vanishes  south  of 
the  auroral-oval,  and  F-layer  propagation  modes  have  their  maximum  frequencies 
reduced.  Also,  the  rays  which  graze  the  auroral -oval  boundary  region  at  F-layer 
heights  experience  the  largest  off -great -circle  deviations  (clockwise  when  viewed 
downward  from  the  top)  from  the  initial  ray  azimuths.  For  initial  ray  azimuths 
lying  south  of  the  auroral-oval  boundary,  and  in  the  trough  region  of  the  critical- 
frequency  function  F^  such  as  shown  in  Figure  6d  for  azimuths  = 50  and  55  degrees, 
and  for  the  transmitter  location  at  (0q,  4>o>  - (32.  0,  1.  8)  degrees,  the  leading  edge 
trace?,  become  simple. 


I.  DISCI  SSION 

At  present,  it  is  thought  likely  that,  in  general,  extensive  ray-trajectory  com- 
putations would  be  required  in  order  to  obtain  the  leading-edge  (minimum  group- 
path)  traces  as  shown  in  the  preceding  section.  This  appears  to  be  true  even  for  a 
relatively  simple  ionosphere  consisting  of  one  single  layer,  if  its  height,  thickness 
and  critical  frequency  all  vary  with  distance. 


8.  Rush,  C.M.,  and  Elkins,  T.J.  (1975)  An  assessment  of  the  magnitude  of  the 
F-region  absorption  on  HF  radio  waves  using  realistic  electron  density  and 
collision  frequency  models,  ITU  Telecommunication  Journal,  1-13. 
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To  infer  the  principal  features  of  an  electron-density  distribution  over  a large 
region  by  means  of  starting  th  a trial  distribution,  and  then  proceeding  by 
iterative  computation  of  lea  liny-edg-  ! races,  a considerable  number  of  successive 
alterations  of  the  initial  trial  distribution  would  need  to  be  made  before  the  com- 
puted leading-edge  traces  car  be  brought  to  acceptable  agreement  with  the  corre- 
sponding traces  in  an  observed  backscatter  ionogram. 

For  this  procedure  of  iterative  computation  of  leading-edge  traces  to  be  gen- 
erally useful,  two  conditions  are  involved.  First,  faster  ray  computation  needs  to 
be  developed.  Now  it  requires  approximately  0.5  hours  or  more  of  Central  Proc- 
essor time  with  the  very  fast  CDC  6000  Computer  in  order  to  obtain  the  leading- 
edge  traces  for  one  initial  ray  azimuth.  Second,  the  procedure  is  workable  pro- 
vided that  when  the  initial  trial  electron  distribution  is  successively  altered  in 
some  monotonic  fashion,  the  resulting  successive  sets  of  leading-edge  traces  are 
also  altered  in  a monotonic  fashion  within  a common  frequency  interval  at  middle- 
range  frequencies  where  the  traces  persist  and  are  relatively  free  from  complicat- 
ing cusps  or  irregularities.  When  this  second  condition  is  met,  multiple-regression 
techniques  can  be  used  to  determine  the  optimal  electron -distribution  being  sought. 
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